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Abstract. Pan-chromatic galaxy surveys are providing tightening constraints on the 
global mass assembly history, and high-resolution imaging of large fields is telling us 
when and where stars were formed. How well are state-of-the-art hierarchical galaxy 
formation models currently doing at reproducing these observations? I present results 
here that suggest that hierarchical models are doing quite well at reproducing the 
global star formation and stellar mass assembly history obtained from galaxies selected 
in the optical and Near IR. However, the same models fail to reproduce two very 
important populations at high redshift: quiescent red spheroids and vigorously star- 
forming, dust-enshrouded starbursts. This mismatch carries important lessons about 
how star formation is triggered and regulated in early galaxies, and may force us to 
consider new ideas about the formation of massive spheroids. 

1 Tracing Galaxy Assembly 

Probing the high-mass end of the galaxy mass function at any epoch provides a 
particularly strong test of theories of galaxy formation. A zoo of apparently mas- 
sive objects at high redshift has been discovered in recent years, in a broad range 
of wavelengths and utilizing a variety of selection techniques. For example, at this 
meeting, we heard about massive galaxies at z 1 in the C0MB0-17/GEMS 
survey (Bell, Rix), VIMOS (Le Fevre), and DEEP (Koo, Newman), near IR se- 
lected galaxies and Extremely Red Objects (EROs) at z ^ 1-2 (Fontana, Daddi, 
Drory, Cimatti), color-selected galaxies in the spectroscopic redshift 'desert' at 
z ~ 1.5-2.5 (Steidel, Chen), and sub-mm and mm selected galaxies at z ~ 2.5 
(Chapman, Bertoldi, Genzel). The challenge posed to us by the multi- wavelength 
universe is to unify this zoo of objects into a coherent picture of galaxy evolution. 

Another exciting observational development is the imaging of relatively large 
fields with the Advanced Camera for Surveys (ACS) on the Hubble Space Tele- 
scope (HST), combined with multi- wavelength ground-based observations. These 
new surveys, such as GOODS (Great Observatories Origins Deep Survey), GEMS 
(Galaxy Evolution from Morphology and SEDs), and COSMOS, will for the 
first time allow us to study the connection between galaxy morphology (through 
structural parameters such as bulge-to-disk ratio) and stellar populations (through 
broad band colors and spectroscopy) for a statistically meaningful sample, over 
a large range of cosmic time. We can then go beyond a global census of star for- 
mation or the build-up of stellar mass, to understand where and how stars were 
formed, and perhaps which processes stimulated or regulated star formation. In 
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this paper, I shall attempt to briefly summarize the status of our theoretical un- 
derstanding of some of these observational results and what we may learn about 
galaxy formation from them. 

2 Do Massive Galaxies at High Redshift Pose a 'Crisis' 
for CDM? 

It has become a familiar story: in the Cold Dark Matter (CDM) model of struc- 
ture formation, small mass objects form first, and larger mass objects form 
hierarchically through mergers and accretion. Therefore, measurements of the 
number density of massive structures at high redshift pose potentially strong 
constraints on this class of models. The extent to which this potential is re- 
alized, however, depends on how massive and how common are the objects at 
any given epoch. Moreover, such arguments also require making a connection 
between an observed population and the dark matter halos in which they are 
expected to reside. For most, if not all, of the populations mentioned above, 
this is far from straightforward. It is a common perception that the detection of 
massive galaxies at high redshift is a serious problem for CDM. It is interesting 
to take a step back for a moment from the complexities of modeling gas physics, 
star formation, feedback, radiative transfer, stellar populations, etc., and to sim- 
ply ask whether the 'concordance' ACDM model produces enough massive dark 
matter halos to plausibly host the objects that have been detected. 

This simple exercise is carried out in Fig. ^ which shows the predicted cu- 
mulative comoving number density of dark matter halos above a given mass, 
from 10^^-10^^ Mq as labeled, as a function of redshift^. These were obtained 
from the Sheth-Tormen modified Press-Schechter model Also shown are 
the estimated comoving number densities of various observed populations, as 
described in the figure caption. It has been claimed that the sub-mm galax- 
ies, EROs and K20 z ^ 2 objects have stellar masses of a few times lO^^M© 
[ 61151231 lU] . Making the plausible assumptions that about fifteen percent of the 
total mass is in the form of baryons, and a few tenths of the baryons are in stars, 
one can easily accommodate the observed numbers of objects within suitably 
massive dark matter halos of ~ 10^'^Mq. The Lyman break galaxies, which are 
more numerous but probably have stellar masses about an order of magnitude 
smaller (a few times IO^^Mq; 27 24 ) are also easily accommodated within halos 
of a few times 10^"'^ M©. Even the z ^ 6 quasars detected by SDSS ^^li which 
must harbor black holes of ~ lO^M© if they are Eddington limited, can be easily 
accommodated within lO^^M© halos (see jSj for more detailed modeling of these 
objects). So far, none of the observed populations poses a fundamental problem 
for concordance ACDM. Had the observed number density of IO^^Mq galaxies 
at z ~ 2 been an order of magnitude higher than the current estimates, theorists 
would have good cause to squirm. As it is, these results suggest that there were 

^ All results presented in this paper assume the 'concordance' ACDM cosmology: fim = 
0.3, Oa = 0.7, Ho = 70 km/s/Mpc, (78=0.9. 
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Fig. 1. Curves show the cumulative number density of dark matter halos more massive 
than 10", 10^^ 10^^, lO", and 10^^ M0(from top to bottom, as labeled), as a function 
of redshift. Points show the estimated comoving number densities of several observed 
populations, as follows. Hexagons: galaxies with stellar masses greater than IO^^Mq 
(lower point) and 2.5 x IO^^Mq (higher point), obtained by integrating the stellar mass 
function of _2, from SDSS-I-2MASS; open square: Extremely Red Objects |22|: ppen 
triangle: K20 galaxies |10|: crosses: sub- mm galaxies 0; six pointed stars: Lyman break 
galaxies I32| . filled square: quasars |13| . For all populations, ACDM predicts that there 
are enough dark matter halos massive enough to plausibly host the observed objects 
(see text). 



fewer massive galaxies in the past, in qualitative agreement with hierarchical 
structure formation. The challenge now lies in understanding how, where and 
why the stars and dust that make these objects visible to us are formed. This 
requires getting our hands a bit more dirty. 



4 R.S. Somerville 




Fig. 2. The global star formation history. The solid curve shows the prediction from a 
semi-analytic model (similar to the 'accelerated quiescent' model of |30p. The dashed 
line shows the model from 0, which was designed to reproduce the counts of bright 
sub-mm galaxies and the Far-IR extragalactic background. Triangles are observational 
estimates based on Ha, squares are based on rest UV, and open diamonds are based 
on mid-IR or sub-mm observations (see |3UI for references). The crosses are new mea- 
surements, based on rest UV, from |16|. Star symbols show measurements based on 
sub-mm observations from 0. All optical/UV measurements are corrected for dust 
extinction and converted to a ACDM cosmology where necessary, as described in |3U|. 

3 The Global Mass Assembly History 

There are several distinct methods of tracing the global history of star forma- 
tion and the build-up of stellar mass in the Universe observationally. One is to 
measure star formation rates by directly counting actively star forming galaxies 
via various tracers (e.g., |2l)l21p . Another is to detect established galaxies at 
different redshifts, estimate their stellar masses from their colors or spectra, and 
sum the masses up to obtain a global average 12 26 17 • Yet a third approach 
is to measure the ages of stellar populations and use galactic 'archeology' to 
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Fig. 3. Global stellar mass assembly history. The solid line shows the total comoving 
density of stellar mass as a function of redshift, from the same semi-analytic model 
shown in Fig. |5] Symbols show observational estimates of the stellar mass density at 
various redshifts from summing the mass in established galaxies. Open circle and many- 
pointed star: local estimates from jS] and Pj. Filled squares: |12| : downward-pointing 
triangles: 'maximum mass' estimates from 12,; crosses: I26| : upward-pointing arrows: 
|17|. The Dickinson et al. and z = points are meant to represent total masses, while 
the Rudnick et al. and Glazebrook et al. points are integrated above a fixed rest-frame 
luminosity or stellar mass, respectively. 

deduce something about star formation at earlier times (e.g. [53]). It is useful 
to compare the results obtained from all three methods, as the attendant errors 
and selection effects are different. 

Fig. 121 shows a compilation of measurements of the star formation history 
of the Universe from z ~ 6 to the present, including the recent results from 
GOODS [121 . The solid curve shows the integrated star formation rate density 
predicted by a fairly standard semi-analytic model (similar to the 'accelerated 
quiescent' model of (SOI)- It was shown in that this model produced good 
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agreement with existing (dust-corrected) observational estimates of the star for- 
mation history to z ~ 4; it is encouraging that the same model agrees almost 
uncannily well with the new measurements out to z ~ 6 from |lf)j . The dashed 
curve is the updated model of P], which I shall return to later. Let us now turn 
to Fig. 13 which shows the total stellar mass density as a function of redshift. 
The solid line is the mass density of all stars from the same semi-analytic model 
shown in Fig. 12 Symbols show observational estimates, as described in the figure 
caption. 

Obviously, in an ideal world, and neglecting recycling from stellar mass loss, 
this curve should be the integral of the curve shown in Fig. 13 For the semi- 
analytic model, this is indeed the case. Given that the semi-analytic model is 
quite a good fit to the observational results shown in Fig. |2 we would obtain 
something quite similar if we simply integrated a curve fit to the observed star 
formation history. However, although if anything the SFR produced by the semi- 
analytic model is a bit on the low side at low redshift (which dominates the 
integral), the total mass produced by the present day is on the high side. At 
z > 1, the observational estimates of the total mass in established systems 
are more than 50 percent lower than the model, or than the integrated star 
formation curve. This could be due to several factors. One possibility is that the 
IMF changes with time (for example, is more top-heavy at early times). Another 
is that the total star formation rate has been overestimated, perhaps because of 
dust corrections, or completeness corrections. A third is that galaxies are missing 
from the stellar mass census (^Hl show that this is likely to be the case at high 
redshift) or that the stellar masses are systematically underestimated. 

This shows that a ACDM-based galaxy formation model has no difficulty 
producing enough early stars or star formation globally. Similar results have 
been found based on hydrodynamic N-body simulations by e.g. Similarly, 
the model predicts that roughly half of the present day stellar mass was in place 
by z ~ 2, and seventy percent by z ~ 1, implying that there should be enough 
old stars (> 8 — 10 Gyr) at z = to satisfy the 'fossil evidence' in present day 
galaxies. 

4 Trouble: Old Red Galaxies and Big Dusty Starbursts 

The new challenge is to go beyond global quantities like those discussed above, 
and to ask what kind of galaxies host those stars and where the star formation 
is occurring. For example, we know that most of the stellar mass in the local 
universe is contained in old, massive, spheroids, while most of the star forma- 
tion is occuring in smaller-mass, disk-dominated objects. The new generation of 
mult i- wavelength surveys is allowing us to determine whether this was true at 
earlier times. The results will pose much more stringent constraints on models. 

Current indications are that the same hierarchical models that do well at re- 
producing global quantities, as discussed above, do much less well at producing 
the correct proportions of different types of galaxies as a function of redshift. In 
particular, while the models seem to produce enough stellar mass globally, at red- 
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shifts of about unity and above they do not produce as many massive/luminous 
red galaxies as are observed. For example, it has been noted in several differ- 
ent studies (e.g. jl4l7l29j ) that semi-analytic models do not produce enough 
'extremely red objects' (EROs), objects with observed-frame optical- infrared 
colors redder than a passively evolving elliptical (e.g. (i? — -fir)vEGA > 5-6), 
which tend to be at redshifts 0.8 < z < 1.5 (see e.g. [12). Semi-analytic models 
underproduce these objects by factors of two to ten, depending on the precise 
color and magnitude cut, and the specific star formation recipes used. Also, al- 
though semi-analytic models can be tuned to give passable agreement with the 
rest frame U — V color-magnitude relation at z = 0, the model fails to produce 
enough luminous/massive, red, bulge-dominated galaxies at z > 0.5-1, based on 
measurements from COMBO-17 and the DEEP redshift survey. 

A possibly related problem is that, when absorption and emission by dust 
are modeled in detail, semi-analytic models also fail to produce enough luminous 
sub- mm galaxies |18lll| . now believed to be predominantly at redshifts z ~ 
2.5 0. This problem is demonstrated in Fig. El the dashed curve shows the 
star formation density needed to account for the sub-mm sources and the Far- 
IR extragalactic background |3I4| . The required star formation at z ~ 1.5-2.5 
greatly exceeds that predicted by the semi-analytic models or hydro simulations. 
One way to solve this problem is to assume that the IMF was very top-heavy at 
early times, or in starbursts (which account for a larger fraction of star formation 
at high redshift), as we heard from Cedric Lacey at this meeting. It is possible, 
though, that these two problems are reflections of the same underlying difficulty: 
perhaps not enough massive objects are being created (in huge dust-enshrouded 
bursts) at z > 2, which is why not enough old, red, quiescent massive objects are 
in place at z ~ 1-2. However, the integrated stellar mass density at low redshift 
poses a strong constraint — it would clearly be exceeded if the Blain et al. curve 
were integrated. This remains a puzzle. 

5 Summary 

To summarize, hierarchical models of galaxy formation (whether based on semi- 
analytic or N-body-|-hydro techniques) are currently able to satisfy observational 
constraints on the global history of star formation and stellar mass assembly from 
optically selected galaxies, from z ~ 6 to the present. However, they fail badly at 
reproducing two very different kinds of massive objects at high redshift: quiescent 
red spheroids, and vigorously star forming, heavily dust-obscured ULIRG-like 
objects detected in the mid-IR and sub-mm. It is possible, though not certain, 
that these problems are connected. One can think of many physical mechanisms, 
neglected in the current modeling, that could quench star formation and produce 
more red objects, but these mechanisms alone could only make galaxies dimmer 
and less massive. What seems to be needed is a mechanism that boosts star 
formation but then rapidly quenches it — in massive systems preferentially. This 
may sound like an old story (monolithic dissipative collapse). What is needed 
is a physical mechanism that produces such a process within the framework of 
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the otherwise largely successful hierarchical structure formation model. There is 
a growing belief that AGN-driven feedback will probably play a major role in 
solving these problems, but many details remain to be worked out. I conclude 
that the multi- wavelength universe is a messy but exciting place to live. 
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